Cotunneling-mediated transport through excited states 
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We present finite bias transport measurements on a few-electron quantum dot. In the Coulomb 
blockade regime, strong signatures of inelastic cotunneling occur which can directly be assigned to 
excited states observed in the non-blockaded regime. In addition, we observe structures related to 
sequential tunneling through the dot, occuring after it has been excited by an inelastic cotunneling 
process. We explain our findings using transport calculations within the real-time Green's function 
approach, including diagrams up to fourth order in the tunneling matrix elements. 

PACS numbers: 73.23.Hk, 73.40.Gk 



In a quantum dot in the Coulomb blockade regime, the 
energy gap related to the charging energy becomes larger 
than k^T and sequential tunneling transport involving 
only dot ground states is exponentially suppressed (see 
e -g- HI)' Transport is dominated by cotunneling Q. 
Elastic cotunneling, prevalent at low bias voltages, in- 
volves virtual tunneling of one electron through the dot 
via a higher-energy state and leaves the dot in the ground 
state. Inelastic processes imply correlated tunneling of 
two electrons, leaving the dot in an excited state with 
energy A above the ground state. Inelastic cotunneling 
sets in once the bias energy eVbi as > A. Recently, Golo- 
vach and Loss have presented a theoretical analysis of the 
interplay between cotunneling and sequential tunneling 
in a double dot system Experimental investigations 
involving cotunneling have been performed on metallic 
El IE EJ an d semiconducting 0, IE 13 systems containing 
a large number of electrons. Signatures of inelastic cotun- 
neling in a transport measurement have been observed in 
investigations on small vertical semiconductor quantum 
dots [loL ITlj , and in single- walled H3 and multiwalled 
|l3j carbon nanotubes, all containing well separated en- 
ergy levels. 

In the following, we first present finite bias transport 
measurements through a quantum dot, showing structure 
outside as well as within the Coulomb blockade regime. 
In the second part, we present a theoretical analysis of 
our results based on transport calculations using the real- 
time Green's function approach. 

The sample (see Fi g. \U a)) was fabricated by surface 
probe lithography [l-l Il5j on a GaAs/Al .3Ga .7As het- 
erostructure, containing a two-dimensional electron gas 
(2DEG) 34 nm below the surface as well as a backgate 
(BG) 1400 nm below the 2DEG. The unstructured 2DEG 
had a mobility of (3.5 ± 0.5) • 10 5 cm 2 /Vs and a density 



of (4.6 ± 0.5) • 10 11 cm~ 2 at a temperature T = 4.2 K and 
a BG voltage Vbg = -0.5 V. 

All measurements were performed in a dilution refrig- 
erator with a base temperature of 80 mK. Negative volt- 
ages were applied to the surrounding gates (see Fig.^a)) 
and to the back gate, to tune the charge on the dot and 
the transparency of its tunnel barriers. The bias voltage 
Vbias was applied symmetrically (with respect to ground) 
across the dot between source (S) and drain (D). The dc 
transport current was measured and numerically differ- 
entiated. An estimated charging energy E c « 1.6 meV 
and a single level spacing A « 0.3 meV were extracted. 

Figure ^ shows measurements of finite bias differen- 
tial conductance dl/dV on a strongly nonlinear scale. 
Figs. IHb)-(e) contain measurements of differential con- 
ductance vs. both gate voltage and bias voltage. Mea- 
surements were performed at different magnetic fields in 
order to vary the wave functions inside the quantum dot 
and their coupling to the reservoirs. Inside the diamond- 
shaped regions, i.e. in the Coulomb blockaded regime, 
we observe horizontal (constant bias) structures. At the 
diamond boundary, the horizontal lines seamlessly join 
some of the most prominent diagonal lines in the non- 
blockaded region. In Figure df), an averaged trace of 
the current vs. bias voltage is presented, showing the po- 
sition of these kinks more precisely. 

For positive bias, e.g. in Fig.nfb), additional structure 
inside the diamond is observed: for bias voltages above 
the well-resolved horizontal threshold line, diagonal lines 
parallel to the diamond edges appear. In our measure- 
ments, this feature remains visible at all magnetic fields 
measured (from B = 0.1 T to B = 0.5 T in steps of 0.1 T, 
see Figs, nfc)-(e) for more examples). The vertical dis- 
tance between the diagonal lines and the diamond edge is 
identical for positively (left of diamond) and negatively 
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FIG. 1: (Color online) (a) AFM micrograph of structure with 
designations of gates: source (S) and drain (D) of the quan- 
tum dot; lateral gates Gl and G2 to control the coupling 
of the dots to the reservoirs; Plunger gates PI and P2 to 
tune the number of electrons on the dot. (b) Finite bias 
measurement of the dot's differential conductance dl /dV at 
B = 0.1 T. (c)-(e) Finite bias transport measurement at 
B = 0.2 T, B = 0.3 T, and B = 0.5 T, respectively. Note 
the different vertical distances between the diamond edge and 
the diagonal structures inside the Coulomb blockaded region, 
(f) Single trace of the bias dependence of the current in the 
blockaded regime, obtained by averaging over the gate volt- 
age range marked in (f) by two vertical dashed lines, (g) 
Color scale used in (b)-(e): To show all features inside and 
outside the Coulomb blockaded region, our data is presented 
in a highly nonlinear color scale. The plotted quantity is cal- 
culated as 1 \/\x\ ■ sign(:r) with x = dl/dV ■ 1 fi. 



sloped lines. When extended towards higher or towards 
negative voltages, most of the diagonal lines apparently 
join prominent lines in the non-blockaded regime. 

A closer examination of the structures reveals a con- 
nection between two energy scales visible inside the 
blockaded region (see Figure|3Ic) for an illustration): an 
extension of each diagonal line intersects the zero bias 
line at a certain point (point A in Fig.^c)). Connecting 
this point to the diamond edge at its intersection with 
the closest horizontal line (point B in Fig. Etc)) yields 
an extension of a diagonal line (of opposite slope) in the 



non-blockaded regime. This remains valid at different 
magnetic fields, where the vertical distance between the 
diagonal line and the diamond edge varies by a factor of 
about 2. 

While we observe the horizontal structures for at least 
12 consecutive Coulomb blockade diamonds, the diagonal 
features are only seen for a maximum of two neighbour- 
ing diamonds up to now. On reducing the transparencies 
of the dot's tunnel barriers, the amplitude of the cotun- 
ncling currents become comparable to the minimum cur- 
rent resolution, and the structures inside the Coulomb 
diamonds gradually disappear. 

We interpret our findings as follows: the horizontal 
lines in the blockaded regime mark the onset of inelas- 
tic cotunneling connected to specific excited states. The 
distance from the zero-bias line corresponds to the single- 
particle level spacing of these states with respect to the 
ground state. At the intersection points at the border of 
the Coulomb diamonds, a direct mapping can be made of 
the excited states that contribute measurably to inelas- 
tic cotunneling and those that open additional transport 
channels in the non-blockaded, finite-bias regime. The 
horizontal line close to zero bias in Fig.^c) suggests the 
presence of a state with low excitation energy contribut- 
ing to inelastic cotunneling. 

The most unconventional features observed are the di- 
agonal lines inside the Coulomb blockaded regions. The 
fact that they have the same slope as the diamond edges 
suggests that they are connected to the alignment of an 
energy level with source (negative slope) or drain (posi- 
tive slope). 

To verify this hypothesis, we have performed trans- 
port calculations within the real-time Green's function 
approach |l6l| including all cotunneling diagrams, i.e. all 
diagrams to fourth order in the tunneling matrix ele- 
ments I R the following, the results for a quantum 
dot with a simple level structure are discussed. 

The simplest quantum dot which shows signatures of 
its excitation spectrum in the Coulomb blockade regime 
is described by an Anderson Hamiltonian Hd having two 
non-degenerate single-particle levels E\ and E 2 and the 
excitation energy A = E2 — E\\ 

H = H B +H R + H T (1) 

H D = E ' c t c i + Uclc lC +c 2 + ea G V G N D (2) 
1=1,2 

The annihilation (creation) operator c ( annihilates 
(creates) an electron of state / in the quantum dot. 
Coulomb interaction is described by the second term of 
(0) leading to an additional interaction energy U when- 
ever the dot is occupied by two electrons. The four pos- 
sible states of the isolated quantum dot are labeled as 
follows: |0, 0) denotes the empty dot, |1,0) the single- 
particle ground state, |1, 1) the excited single-particle 
state, and |2,0) the two-particle state. Here we assume 
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that the applied gate voltage Vq leads to a constant elec- 
trostatic potential described by the third term of the 
quantum dot Hamiltonian. In this term, N^> = cfci 
is the number operator for the dot electrons and qq the 
electrostatic lever arm of the gate electrode. The cou- 
pling of the quantum dot to two reservoirs is described by 
the reservoir (Hr) and the tunneling Hamiltonian (Ht)- 
They are of the conventional form (see e.g. [J3 > E3 ) with 
the reservoir electrons being treated as non-interacting 
except for an overall self-consistent potential owing 
to the high density of states in source and drain con- 
tacts. For simplicity, we assume for the following the 
absolute value of the complex tunneling matrix elements 
to be independent of all quantum numbers with a com- 
plex phase which is random with respect to the direction 
of the reservoir electrons wave vector. 

In order to calculate the non-equilibrium trans- 
port properties for finite transport voltages Vsd, 
we use the real-time transport theory developed by 
Schoeller et al. [lfij . Following the steps of this theory 
one can trace out the reservoir degrees of freedom and 
derive a formally correct equation of motion for the re- 
duced density matrix of the quantum dot system which 
under steady state conditions transforms into 

l(E s -E' s )Pf s , = J2PT lS - f dt'X ss , SlS ,(0,t'). (3) 

J -°° 

Here P ss ' denotes a matrix element of the reduced dot 
density matrix with the (few)-particle states \s) and \s') 
of the isolated quantum dot. The kernel S ss / SlS ^ repre- 
sents a generalized transition rate involving the relevant 
tunneling processes. Within the same formalism one can 
also calculate the tunneling current expectation value for 
the steady state 

(Z?) =-e£ P^, f dt'Z: s8iS ,(0,t'). (4) 

/ J — oo 

In the sequential tunneling approximation, a "tunnel- 
ing in" process is only possible if a reservoir electron 
matches the energy required to charge the quantum dot 
by a further electron. Generally, this energy is given by 
E s — E s i for a transition between the state \s) (where 
\s) = \N,j) is the jth TV-particle state) and the (N + 1)- 
particle state |s'} and is in the following called trans- 
port channel and denoted by /x(s;s'). For the quan- 
tum dot described by @, four transport channels ex- 
ist. In the Coulomb blockade regime, where transport 
in lowest order is exponentially supressed, the electro- 
chemical potentials of both reservoirs are in between the 
two transport channels associated with two groundstates: 
/i(2,0;l,0) > fj. r > 0; 0, 0). Going one step fur- 
ther and calculating the kernel of {3} in fourth order, 
already 64 qualitativly different terms occur, describing 
so-called cotunneling processes in which two electrons 




ea G V G /kT dl/dV SD (a.u.) 



FIG. 2: (Color online) (a) Calculated charging diagram within 
the Coulomb blockade regime including the cotunneling con- 
tribution. The peak conductances due to resonant sequential 
tunneling (also responsible for the increasing background in 
(b)) have been "cut" from the color scale (outer white re- 
gions). Solid, dotted and dashed lines correspond to traces in 
(b). (b) Cotunneling signatures within the Coulomb block- 
ade regime. Shown is the differential conductance versus the 
applied transport voltage for three different gate voltages. 
The arrow marks the onset of inelastic cotunneling for the 
dashed trace. The parameters used for this calculation are: 
T = 0.1/sbT, S = 14.67fc B T, and U = 52fc B T. Bottom: color 
scale for (a). 

participate coherently in a tunneling process. In the fol- 
lowing we study the resulting differential conductance 
including consistently all cotunneling contributions |2lj . 

In Fig. Hfb) the differential conductance as a function 
of the applied source-drain voltage Vsd is shown within 
the Coulomb blockade regime, i.e., the electro-chemical 
potentials for Vsd = are energetically in between the 
two ground state channels /i(2,0;l,0) and //(l, 0; 0, 0). 
For small voltages, all three traces start with the same 
value and at least the traces for eaoVc / 'k^T = 29.3 and 
eaoVQ/k^T = 33.3 stay constant for small transport 
voltages. This constant and finite differential conduc- 
tance can be attributed to elastic cotunneling by virtual 
tunneling through either the vacuum or the two-particle 
state. Additionally, for all three traces a peak is found 
which shifts linearly to higher source-drain voltages with 
increasing gate voltage. In contrast to the traces at lower 
gate voltages, the differential conductance of the highest 
gate voltage (eaQVc/fceT = 33.3) shows an additional 
step (see arrow in Fig.^b)) emerging at the source-drain 
voltage eV$T>/k,B,T « 15, which corresponds to the exci- 
tation energy S. This step is strongly smeared due to 
temperature and especially due to the overlap with the 
peak occurring at higher voltages (also responsible for 
the strong increase towards higher bias voltages). 

In Fig.[3fa), the relative position of the transport chan- 
nels with respect to the electrochemical potentials in the 
reservoirs is shown at the parameters where the step oc- 
curs (eVsD = 8 and eaaVc/k^T = 33.3). In this situa- 
tion, inelastic cotunneling becomes possible in which the 
quantum dot becomes excited during the tunneling pro- 
cess. This additional (to the elastic process) cotunneling 
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FIG. 3: (Color online) (a/b) Relative position of transport 
channels and contact electrochemical potentials for different 
bias voltages, (c) Schematic picture of the various tunneling 
regimes within the Coulomb blockade. A and B label the 
marked intersection points, (a) and (b) mark the gate voltage 
values of the corresponding diagrams. E is for elastic and I 
for inelastic cotunneling, S for sequential tunneling. 



process leads to a step in the differential conductance 
[20| which corresponds to the experimentally observed 
steps showing as horizontal lines in Fig. 1. Eventually, 
by further increasing the source-drain voltage eVsD > 8, 
the electrochemical potential of the drain reservoir be- 
comes resonant with the transport channel /i(l,l;0,0) 
(sketched in Fig. EJb)). Other than in the sequential 
tunneling approximation, where the |1,1) state cannot 
be occupied due to the Coulomb blockade effect, inelas- 
tic cotunneling allows to occupy this excited state and 
the resonant channel leads to a peak in the differential 
conductance. Due to the smaller cotunneling rate, the 
peak is lower as compared to the corresponding peak be- 
yond the Coulomb blockade regime. 

For lower gate voltages, the peak moves to lower 
source-drain voltages and eventually merges with the step 
at eVsD = 8. For even lower gate voltages, the channel 
fi(l, 1; 0, 0) is already within the transport window at the 
source-drain voltage eVsD = 8 needed to allow for the in- 
elastic cotunneling process. 

Combining all these processes, various tunneling 
regimes within the Coulomb blockade can be identified 
as sketched in Fig.^c). For eT^D < 8, transport is dom- 
inated by elastic cotunneling, leading to a constant offset 
of the differential conductance. For gate-voltages in the 
vicinity of the Coulomb blockade center and eVsD > 8, 
a regime where elastic and inelastic cotunneling occur is 
found. In the remaining outer regime, sequential tunnel- 
ing through the excited single-particle state is also pos- 
sible. At the border of this regime, a peak occurs in 
the differential conductance. All described features are 
also found in the calculated charging diagram including 
cotunneling (shown in Fig. E{a)). 

While it is not entirely clear to us why the induced se- 
quential tunneling contributions have not been observed 
before, we can identify a few requirements: First of all, 
the charging energy has to be large enough, E c > 28. 
This can be directly seen from Fig. EIc), and it explains 



why the effect has not been observed e.g. in [lj}. I n 
addition, a sufficient level spacing 8 » k^T is required 
so that the effect is not smeared out by temperature. 
We speculate that the strongly asymmetric tunnel bar- 
rier configuration used in the present experiment may 
also enhance the visibility of these features. 

Elastic cotunneling has earlier been identified as a pos- 
sible source of uncertainty in the operation of single- 
electron devices (see e.g. 0). Our results show that 
the inelastic contributions can become more prominent, 
especially if the induced sequential tunneling is taken 
into account. It follows that in an application relying 
on Coulomb blockade in quantum dots, e.g. in quantum 
information processing, the bias must be kept small in 
comparison to the lowest excitation energy. 
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